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Abstract 

A technique was developed to calculate the stress 
intensity factor for multiple interacting cracks. The 
analysis was verified through comparison with ac- 
cepted methods of calculating stress intensity fac- 
tors. The technique was incorporated into a fatigue 
crack growth prediction model and used to predict 
the fatigue crack growth life for multiple-site damage 
(MSD). The analysis was verified through compari- 
son with experiments conducted on uniaxially loaded 
flat panels with multiple cracks. Configurations with 
nearly equal and unequal crack distributions were ex- 
amined. The fatigue crack growth predictions agreed 
within 20 percent of the experimental lives for all 
crack configurations considered. 

Introduction 

The continued structural integrity of the aging 
commercial transport aircraft fleet is of great con- 
cern to the aerospace community. The long service 
life of these aircraft increases the possibility of a re- 
duction or loss of structural integrity due to fatigue 
cracking. Multiple-site damage (MSD), in particular, 
1 of fatigue cracking that presents special 
o the aircraft maintenance operator as well 
,ructural analyst. MSD refers to the occur- 
3veral cracks close enough together to in- 
;h other and to affect the overall structural 
The critical crack size for individual cracks 
latively small, making their detection with 
ondestructive examination (NDE) methods 
The mutual influence of the adjacent cracks 
the complexity of predicting fatigue crack 
3havior. 

near elastic fracture mechanics (LEFM) ap- 
i predicting fatigue crack growth relates the 
snsity factor to the crack growth rate. Stress 
factors can be calculated with a variety of 
,1 and numerical techniques. Closed form so- 
snerally exist for a limited number of simple 
:h as infinite plate problems (ref. 1). Nu- 
;echniques such as finite element or bound- 
lent analyses are often used for more com- 
:k configurations (ref. 2). However, the large 
of crack tips potentially involved in analyzing 
1 problem would require extensive computer 
s using finite element analyses. Boundary 
dmiques (ref. 3) require less modeling effort, 
only model one crack. Approximate solu- 
MSD stress intensity factors were developed 
oining known solutions (refs. 4 and 5). 

rerical techniques could be used to calculate 
ntensity factors for the entire fatigue crack 


growth life. However, the scale of modeling required 
to capture the effect of both the local structural 
details and the large-scale structural details would 
be prohibitive in a life analysis where the stress 
intensity factors are calculated many times as the 
cracks grow. A more practical approach would be to 
consider what influences the crack growth behavior 
at different stages of life. The damage tolerant life of 
a structural component susceptible to MSD cracking 
can be divided into three regions of crack growth: 
local, MSD, and post-linkup. The local crack growth 
region consists of the initial development of cracks at 
a local structural detail, such as a rivet hole. The 
cracks are strongly influenced by the local structural 
details, but are too small to have a large influence on 
adjacent cracks. The cracks in the MSD region are 
strongly influenced by the adjacent cracks. The crack 
in the post-linkup region behaves as a single long 
crack influenced by the large-scale structural details. 

The technique described in this paper examines 
the crack growth life for the local and MSD regions 
of crack growth. The post-linkup region was not 
considered and would require large-scale analyses of 
structural details, such as stringers, cutouts, and 
riveted connections, using finite element shell codes. 
An alternating indirect boundary element (AIBE) 
technique was developed herein to calculate the stress 
intensity factor for multiple interacting cracks. A 
hole correction factor was added to the AIBE solution 
to reduce the detail of modeling required. This 
technique was sufficiently fast and accurate to be 
incorporated into a fatigue crack growth prediction 
code. 

This paper describes the AIBE technique of de- 
termining stress intensity factors. The technique was 
verified by comparisons of the calculated stress in- 
tensity factors with those determined from accepted 
solutions, finite element methods, and boundary el- 
ement methods. Damage tolerant life predictions 
were made with the AIBE technique in conjunc- 
tion with a hole correction factor. The life predic- 
tions were compared with MSD tests conducted on 
flat 2024-T3 aluminum alloy specimens with multiple 
cracks propagating from a line of open holes. 

Symbols 

a half crack length 

BFM boundary force method 

C constant 

d half centerline distance between 

holes 

dj half element length 



K 


S 


stress intensity factor 
Np number of cycles, predicted 

Nf number of cycles, test 

R stress ratio 

S applied stress 

s integration variable 

w half specimen width 

a constraint 

(j 0 flow stress 

a n normal stress 

r shear stress 

Stress Intensity Factor Calculations 

The use of linear elastic fracture mechanics to pre- 
dict fatigue crack growth behavior requires the de- 
termination of the stress intensity factor. The stress 
intensity factor is influenced by geometric consider- 
ations such as crack length, other adjacent cracks, 
free boundaries, and holes. The various geometric 
considerations will affect the stress intensity factor 
differently as the crack grows. Consider, for exam- 
ple, cracks propagating out of two open holes in an 
infinite body, as shown in figure 1. The influence of 
:iie hole is large when the ratio of half crack length 
to hole radius a/r is small and then decreases as a/r 
increases, while the influence of the adjacent crack 
steadily increases, as shown in figure 2. 

The growth of multiple cracks in a structure 
can be divided into three regions, each of which is 
governed by different geometric considerations. The 
stress intensity factor of cracks in the local crack 
growth region is strongly influenced by the local 
structural details. The cracks in this region are 
small, as would be the influence of any small adjacent 
cracks. The stress intensity factors of cracks in the 
MSD crack growth region are strongly influenced by 
the adjacent cracks, and the effects of local structural 
details diminish. The post- linkup region is influenced 
strongly by large-scale structural details. 

An alternating indirect boundary element (AIBE) 
method was developed to calculate the stress inten- 
sity factor for multiple interacting cracks. A correc- 
tion factor was applied to the stress intensity factor 
to account for the local structural details. Consid- 
eration of the post-linkup region was not required 
for the structural configuration examined. The fol- 
lowing sections describe the AIBE technique and the 
hole correction factor. 



Figure 1. Cracks propagating from two sides of two holes in 
an infinite plate subjected to uniaxial loading. 



Figure 2. The influence of a hole and an adjacent crack on the 
stress intensity factor for the cracks illustrated in figure 1 
(r/d = 0.15). 


Alternating Indirect Boundary Element 
(AIBE) Method 


The stress intensity factors for multiple interact- 
ing cracks were calculated with an alternating in- 
direct boundary element (AIBE) technique. The 
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AIBE method is based on the principle of superpo- 
sition of stresses associated with cracks and free (or 
loaded) boundaries. To illustrate the superposition 
of stresses, consider a single crack in a semi-infinite 
sheet, as shown in figure 3. (An MSD configura- 
tion would be obtained by adding additional cracks 
and free boundaries.) Figure 4 illustrates how su- 
perposition was used to determine the stress func- 
tions for a crack in a semi-infinite sheet. The cracked 
infinite sheet with a uniform crack-face stress S is 
shown in figure 4(b). The uncracked infinite sheet 
subjected to a uniform remote stress S is shown in 
figure 4(c). Figure 4(d) shows the uncracked sheet 
with the nonuniform stresses cr b x (y) and r^ y (y) act- 
ing on a line that corresponds to the free boundary 
and the nonuniform stresses (T%{x) and r£ y (x) act- 
ing along y — 0, |x| < a. The crack-face loading 
(fig. 4(e)) is superimposed on the uncracked loading 
(fig. 4(d)) to produce the stress-free crack face shown 
in figure 4(a). 



Figure 3. Semi-infinite body with crack. 



(a) Cracked semi-infinite sheet 
with remote stress S. 


(b) Cracked infinite sheet with 
uniform crack-face stress. 


(c) Uncracked infinite sheet with 
uniform remote stress. 




(d) Uncracked infinite sheet with nonuniform (e) Infinite sheet with crack-face loading, 

stresses acting along free boundary and along 
V = 0, \x\ < a. 


Figure 4. Principle of superposition applied to a crack in a semi-infinite plate. 


The formulation requires analytical expressions 
for the stresses due to the loadings shown in fig- 
ures 4(d) and 4(e). Analytical expressions developed 
by Timoshenko and Goodier (ref. 6), Westergaard 
(refs. 7), and Tada and Irwin (ref. 8) were used to 
formulate the necessary expressions for the stress dis- 
tribution. 


P x 

7y 2 tt x 2 4- y 2 

_ -P y 

xy 2n x 2 + y 2 


(1 -v) -2(l + i/) 


x 2 + y 2 


(1 - v) + 2 (1 + v) 


x 2 + y 2 


( 2 ) 

( 3 ) 


The stresses due to a force Q acting in the y-direction 
are given by 


Uncracked infinite sheet with nonuniform 

line loading. The analytical expressions defined by 
Timoshenko and Goodier (ref. 6) for the stresses at 
any point {x,y) due to a concentrated force on an 
infinite body were used to formulate the equations 
for the stresses in an uncracked infinite sheet with 
a nonuniform line loading (fig. 4(d)). In the Tim- 
oshenko and Goodier expressions, the concentrated 
force is applied at the origin of the local coordinate 
system (x ,y), as illustrated in figure 5. The stresses 
due to a force P acting in the x-direction are given 
by: 
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( 5 ) 

( 6 ) 


where u is Poisson’s ratio. The stresses due to normal 
and shear stress distributions were obtained by divid- 
ing the boundary into discrete elements and integrat- 
ing equations (l)-(6) over each element length (2 dj), 
as illustrated in figure 6. The stresses due to a normal 
stress distribution cr n (s) are given by 
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where a n (s) is assumed to have a quadratic form as 


&n (s) = Co + C\S + C2i 


The stresses due to a shear stress distribution r(s) are given by 
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Figure 6. A line element with a quadratic normal stress 
distribution for an infinite body (note that a similar 
procedure is used for the shear stresses). 


where r(s) is assumed to have a quadratic form as 

t (s) = Do + D\s + D2S 2 (14) 

where Ci and Di ( i = 0, 1, 2) are known for each ele- 
ment and will be discussed later. A coordinate trans- 
formation is performed on the calculated stresses 
to be consistent with the global coordinate system 
(X,Y). 

Cracked infinite sheet with nonuniform 
crack-face loading. From Westergaard (ref. 7), the 


stresses in an infinite sheet with crack- face loadings 
can be written in terms of two stress functions, Zj 
and Zjj , as follows: 

cr x — Re Zj — y lm + 2 lm Zji + y Re ( 15 ) 

dz ciz 

<T y = Re Zj + y lm - — y Re — j— (16) 

dz dz 


Txy = -y Re + Re Z n - y lm (17) 


The Westergaard stress functions, Zj and Zjj, are 
defined for point forces P and Q acting on a crack of 
length 2a, as illustrated in figure 7. 



Figure 7. Concentrated forces P and Q acting on a crack in 
an infinite body. 


Zi (*) = 

P 

Va 2 - b 2 

(18) 

7 r 

(z — b) V z 2 — a 2 

Z n (z) = 

Q 

Va 2 - b 2 

(19) 

7 r 

(z — b)V z 2 — a 2 


where z = x + iy. 

The stress functions for crack-face stress distri- 
butions were obtained by dividing the crack into 
discrete elements and integrating equations (18) 
and (19) over the length of each crack element (2 dj), 
as illustrated in figure 8. The resulting stress 
functions for each element are given by 




Figure 8. A crack element with a quadratic normal stress distribution for a crack in an infinite body (note that a similar procedure 
is used for the shear stresses). 



Equations (19) and (20) could theoretically be in- 
tegrated exactly, but numerical integration proved 
to be computationally efficient. A coordinate trans- 
formation was performed on the calculated stresses 
to be consistent with the global coordinate system 

pen 


Stress intensity factor equations. The 
mode I and mode II stress intensity factors due to 
concentrated loads P and Q on a crack in an infinite 
body (refs. 7 and 8) are given by 




The stress intensity factors for a crack face stress 
distribution are determined by integrating equa- 


tions (22) and (23) over the length of each crack 
element (2 df) measured with respect to the dis- 
tance from the element centroid to the center of the 



(24) 


The total stress intensity factor solutions are ob- 
tained by adding the contributions from each crack 
element. 


MSD configurations. A two-dimensional 
cracked structure can be modeled with the alternat- 
ing indirect boundary element (AIDE) method by de- 
scribing the cracks and external boundaries with a 
series of line segments or elements. The boundaries 
must be continuous and the cracks contained within 
the boundaries and not intersecting other cracks. 
The quadratic elemental stress distribution was used 
to reduce the number of elements (i.e., the number 
of degrees of freedom) required to describe a crack 
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problem. The quadratic stress distributions require 
that the stresses be defined at three nodal points 
in each element, s = —0.6 dj, 0, and 0 .6dj. The 
solution is obtained by first assuming initial values 
of the six stress coefficients (eqs. (10) and (14)) of 
each element. The initial value is not crucial, a 
value of unity for each coefficient is sufficient for 


convergence, but using the results from a previous 
solution (when available) will considerably reduce 
the number of iterations. The stresses for each 
element (three normal and three shear stresses) are 
determined with the current stress coefficients and 
assembled into a vector of length 6n, where n is the 
total number of elements, as 





r(0.6dj_i) 


D^ 1 + DI 1 (0.6di_ x ) + D{~ 1 (O.Gd^i) 2 

On (—0.6 di) 


q + C[ (-0 ,6di) + C\ (-ORdi) 2 

r(-O.Qdi) 


D'q + D\ (— 0 Mi) + D\ (-0.6d, ; ) 2 

°n{ 0) 

. = < 

ci o 

r( 0) 

n?' 

o n (0.6 di) 


q + C[ (0.6di) + C\ (OM) 2 

t (0.6di) 


Dj + D{ (0.6 di) + Dj (0.6di) 2 

o n (-0.6dj + i) 


q +1 + c\ +l (—0.6 di +1 ) + q +1 (-0.6 d i+1 f 


• / 


An influence matrix [t] is formed using equa- 
tions (7) — (17) to determine the effect of the stresses 
{o- old } of element i on the stresses {cr new } at ele- 
ment j. The influence matrix is a fully populated, 
nonsymmetric 6n x 6n matrix. The new elemen- 
tal stresses are defined such that the imposed stress 
constraints { CT im P° s edj. are satisfied, as given by 

{a 1 *™} = ^imposed J _ ^ j^old j ( 2 7) 

The imposed stress constraints {<r imposed } are ei- 
ther zero for stress-free boundaries and cracks or 
the applied stresses. Once the stress distribu- 
tions {cr new } are calculated, the stress intensity 
factors are calculated with equations (24) and (25). 

The process is repeated until the stress inten- 
sity factor solution converges, typically defined as 
changes of less than 0.5 percent. The solution could 
be obtained directly by inverting [f], as given by 

{<7 new } - [r 1 ] |cr imposed | (28) 

However, the iteration method was more efficient for 
use in a life prediction code, which requires repeated 
stress intensity factor solutions for a large number of 
crack lengths. Rapid convergence of the solution is 
obtained if the stress coefficients from the previous 


crack lengths are used as the initial values in the 
current calculation. 

Local Correction Factor 

The stress intensity factor for cracks in the local 
growth region is strongly influenced by local struc- 
tural details. To describe these details explicitly 
would add considerably to the size of the required 
MSD model. Alternatively, one structural detail 
could be modeled accurately, and a simple relation- 
ship describing its influence on the stress intensity 
factor could be developed. Such a relationship could 
be incorporated into the MSD stress intensity factor 
calculations. 

One structural detail that has been extensively 
analyzed is an open hole. The stress intensity factor 
for cracks propagating from opposite sides of an open 
hole in an infinite plate can be written as 

K = S^F h 0 (29) 


= Jl-~ 1.0 + 0.358^ + 1.425 (-Y 


1.578 ( -) 3 + 2.156 


(30' 





Figure 9. Illustration of use of the hole correction factor F h ( p) for cracks emanating from circular holes. 


where r is the hole radius, S is the applied stress, 
a is the half crack length, and Ff, (®) is the hole 
correction factor (ref. 9) . The stress intensity factors 
for a series of cracks emanating from a row of holes 
would be the products of the hole correction factor 
for each hole and the stress intensity factor for each 
crack, as illustrated in figure 9. 

Stress Intensity Factor Verification 

The AIBE technique was verified by comparing 
the stress intensity factors calculated by AIBE with 
those from an accepted solution, a boundary force 
analysis, and finite element analyses. The crack 
configurations used both a single center crack and 
multiple cracks in a finite body. The use of the local 
correction factor was examined by considering cracks 
propagating from two holes. 

AIBE Stress Intensity Factor Verification 

A center crack tension (CCT) panel subjected to 
uniaxial loading perpendicular to the plane of the 
crack was analyzed. The crack was divided into 
12 elements, with smaller element lengths near the 
crack tip, as shown in figure 10. The boundary 
mesh consisted of 72 elements, with smaller element 
lengths near the corners and the edges closest to 
the crack tip. (The symmetry of the configuration 
was not used to reduce the model.) The variation 
in element sizes was necessary to account for the 
steep gradients in the required balancing stresses that 
occur at points of discontinuity in the boundary and 
are not necessarily associated with steep gradients 
in physical stresses. An accepted solution for the 
stress intensity factor for a CCT specimen is given 
by Fedderson in reference 10 as 


K = sjnase c (~) (31) 

The stress intensity factor results from the AIBE 
analysis agreed with the accepted solution within 
1 percent for crack length- to-- width ratios (a/w) 
below 0.8, as shown in figure 11. 



Figure 10. AIBE mesh used in the analysis of a center crack 

tension specimen, h/w = 2. 

A centered, inclined crack subjected to uniaxial 
loading was also analyzed. The crack length was 
constant with a/w = 1/3, and the orientation of the 
crack with respect to the direction of loading was 
varied from 9 = 0° to 85° ( 9 = 0° being perpendic- 
ular to the direction of loading). The mode I and 
mode II stress intensity factors were calculated with 
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Figure 11. Normalized stress intensity factor results for a 

CCT specimen, h/w = 2. 

the AIBE code and the BFM boundary force code 
(ref. 3). The AIBE boundary mesh was the same as 
that used in the analysis of the CCT configuration, 
and the crack mesh consisted of 12 elements oriented 
at an angle 9. The BFM mesh contained roughly the 
same number of elements. The two methods agreed 
within 0.7 percent for the mode I stress intensity fac- 
tor and 0.3 percent for the mode II stress intensity 
factor, as shown in figure 12. 

The next configuration analyzed was that of two 
collinear cracks in a finite plate. The cracks were 
oriented perpendicular to the direction of loading 
and located along the centerline of the specimen, as 
shown in figure 13. The length of both cracks was 
held constant ( a/w = 0.0625) and the distance be- 
tween them (2d) was varied. The stress intensity 
factors were calculated with the AIBE code and the 
FRANC finite element code (ref. 11). The AIBE 
boundary mesh was the same as used in the earlier 
examples, and each crack mesh consisted of 12 ele- 
ments. Because of symmetry, the finite element mesh 
modeled one quarter of the specimen and consisted 
of 475 6-noded triangular elements. The stress inten- 
sity factor results of the two methods agreed within 
1 percent, as shown in figure 14. 


Stress Intensity Factor Hole Correction 
Verification 

The use of a correction factor for the stress in- 
tensity factor was examined by comparing the stress 
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Figure 12. Normalized stress intensity factor for a single 
centered crack in a finite body, with the crack orientated 
at an angle 9 (6 = 0° being perpendicular to the direction 
of loading) h/w — 2; a/w = 1/3. 
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Figure 13. Illustration of multiple crack configuration used 
in the stress intensity factor verification, h/w = 2; 
a/w = 0.0625. 
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by the hole correction factor for cracks propagating 
from an open hole (eq. (29)). The stress intensity 
factor for either of the two inner crack tips calcu- 
lated by the AIBE analysis with and without the hole 
correction factor is compared with results from the 
FRANC finite element code (ref. 11), as shown in fig- 
ure 16. Because of symmetry, the finite element mesh 
modeled one quarter of the specimen and consisted 
of 444 6-noded triangular and 8-noded quadrilateral 
elements. The stress intensity factor results from the 
AIBE method with the hole correction factor agreed 
within 1 percent with the finite element results, as 
shown in figure 16. 



Figure 16. Stress intensity factor for the inner two crack tips 
shown in figure 15. 

Experimental Procedure 

MSD fatigue crack growth tests were conducted 
on five 304-mm-wide 2.29-mm-thick 2024-T3 alu- 
minum alloy specimens subjected to remote uniform 
stress, as shown in figure 17. The specimens con- 
tained 10 holes, which were 3.8 mm in diameter and 
spaced 25.4 mm apart on centers. Cracks were de- 


Figure 15. Illustration of the multiple open hole with cracks veloped by cutting small notches (0.7 mm long by 

configuration used in the stress intensity factor verification. 0.2 mm high) and cycling at the intended fatigue 

r /d = 0.15; r/w = 0.0125. stress until distinct fatigue cracks were propagat- 


intensity factor calculated from the AIBE technique 
with the hole correction factor with results from a 
finite element analysis. The configuration considered 
was that of cracks propagating from two holes in a 
finite plate, as shown in figure 15. The AIBE pro- 
gram calculated the stress intensity factor for the two 
cracks in a finite plate. This solution was multiplied 


ing from both sides of each hole. The experiments 
were conducted under constant amplitude loading 
(71 MPa) at a stress ratio of R = 0.0 and a fre- 
quency of 5 Hz. In three specimens (A6-02, A6-04, 
and A6-05) the center two holes were connected by a 
saw cut after precracking, as illustrated in figure 18, 
to simulate linkup of two adjacent cracks. The initial 
crack lengths after precracking are given in table I. 
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Figure 17. Schematic of the MSD cracking configuration. 



Figure 18. Illustration of saw cuts used to cause premature 
linkup of two adjacent cracks. 


Life Prediction Procedure 

Fatigue crack growth predictions require a de- 
scription of the stress intensity factor and the 
baseline crack growth rate behavior. The stress in- 
tensity factors used in the predictions were obtained 
from the AIBE analysis with the correction factor for 
an open hole in an infinite plate for the configuration 
shown in figure 17. 


Table I. Measured Initial Half Crack Length for the 
Constant Amplitude Fatigue Tests 


Hole 

number 

Crack 

tip 

Crack length, mm 

A6-01 

A6-02 

A6-03 

A6-04 

A6-0-5 

1 

O.L 

2.90 

3.61 

2.84 

2.82 

3.48 

1 

a R 

2.87 

3.02 

2.62 

3.23 

2.95 

2 

dL 

3.25 

4.52 

3.00 

3.63 

5.08 

2 

dR 

3.25 

4.19 

2.31 

4.34 

4.09 

3 

a L 

3.63 

3.30 

3.99 

4.75 

3.23 

3 

dR 

3.63 

3.45 

3.96 

5.18 

2.84 

4 

a-L 

2.62 

3.48 

3.51 

5.16 

3.15 

4 

dR 

2.62 

3.76 

2.90 

5.13 

3.05 

5 

dL 

3.15 

16.0 

3.76 

17.0 

16.2 

5 

dR 

3.15 

* 

2.34 

* 

* 

6 

dL 

2.57 

* 

2.41 

* 

* 

6 

dR 

2.57 

16.0 

3.33 

17.0 

16.2 

7 

dL 

2.39 

3.15 

4.62 

3.51 

3.25 

7 

dR 

2.39 

2.74 

4.06 

3.53 

3.45 

8 

dL 

2.46 

2.69 

2.77 

3.96 

3.28 

8 

dR 

2.44 

2.82 

3.48 

4.04 

3.20 

9 

dL 

3.30 

3.45 

3.84 

3.89 

3.73 

9 

dR 

3.30 

3.76 

5.08 

3.78 

3.94 

10 

dL 

3.00 

2.79 

2.97 

2.95 

3.10 

10 

dR 

3.00 

3.30 

3.51 

3.43 

3.81 


‘Linkup between adjacent cracks. 
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The baseline fatigue crack growth rate behavior 
for 2024-T3 aluminum used in the crack growth 
predictions was obtained from references 12 and 13 
from constant amplitude tests on center crack tension 
specimens. The effects of fatigue crack closure were 
taken into account by calculating the crack opening 
stress S 0 from an equation (ref. 14) based on a 
two-dimensional closure model (ref. 15), as given by 

= A q + A L R + A 2 R 2 + A 3 R 3 (R > 0) (32) 

‘-’max 

-—-^Aq + AxR (R< 0) (33) 

‘ J max 


Table II. Effective Stress Intensity Factor Range Against 
Crack Growth Rate for 2024-T3 (refs. 12 and 13) 


AK eg , 

da/dN , 

MPa-m 1 / 2 

m/cycle 

1.42 

3.56 x 10~ 10 

2.41 

3.00 x 10~ 9 

3.29 

6.10 x 10~ 9 

6.18 

1.00 x 10~ 7 

10.95 

4.32 x 10“ 7 

27.38 

1.78 x 10~ 5 

49.29 

2.54 x 10“ 4 


(34) 

(35) 

(36) 

(37) 

where 

S'rnax maximum applied stress 

R stress ratio (minimum stress/maximum 

stress) 

cr 0 flow stress (average between uniaxial 

yield stress and uniaxial ultimate tensile 
strength, 410 MPa for 2024-T3) 




(0.825 - 0.34a + 0.05a 2 ) cos ^^ --)] ' 


dj. = (0.415 - 0.071a) 


A‘2 = 1 -- do — Ai - d 3 
d 3 = 2 d 0 + A x - 1 


The constraint factor a simulates the three- 
dimensional constraint. A value of a — 1.73 was 
found to be appropriate for 2024-T3 aluminum alloy 

(B = 2.29 mm) (ref. 15). 

The crack opening stress, given by equation (32), 
was used to calculate an effective stress intensity 
factor range for the baseline crack growth rate data. 
The effective stress intensity factor range is given by 


AA' e ff = 


A K 
l-R 



(38) 


Crack growth rate da/ dN and effective stress inten- 
sity factor range A K e g data were determined for use 
in the life predictions, as given by table II. 

The life prediction started with the initial dis- 
tribution of the 20 crack tips. The stress intensity 
factor and AK e fi were determined for each crack tip, 
and the crack growth rate da/dN obtained from the 
baseline behavior (table II). The crack tip with the 
largest crack growth rate ( da/dN) max was grown a 


fixed increment (A a = 0.12 mm), and the required 
number of applied load cycles AN was determined 
as follows: 


AN = 


A a 

da/ dN max 


(39) 


The crack growth increments for the remaining crack 
tips were determined from the number of applied load 
cycles as 


A ai = 


da 

dNi 


AN 


(40) 


where i is the crack tip number. The process was 
repeated until the net section stress was equal to 
the yield stress. The increments of applied load 
cycles AN were added to obtain the predicted cycles 
to failure N p . 


Life Prediction Results 

The MSD fatigue tests were cycled to failure 
with optical crack length readings made periodically 
throughout the tests. The average initial half crack 
length and the observed cycles to failure are given 
in table III. The average half crack length was the 
sum of the tip-to-tip crack lengths divided by 20, 
regardless of whether any linkup of the cracks had 
occurred. 


Table III. Summary of Constant Amplitude Fatigue Tests 


Specimen 

ID 

Average 
initial 
half crack 
length, 
mm 

N t 

Cycles 

to 

failure 

A6-01 

2.926 

58 740 

A6-02 

*4.303 

26430 

A6-03 

3.368 

46210 

A 6- 04 

*4.864 

15 600 

A6-05 

*4.400 

25 920 


‘Holes 5 and 6 were connected by a saw cut. 
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As described above, life predictions were made 
for each specimen by using the baseline crack growth 
rate behavior given in table II, AA e fj- calculated 
with equation (38), and the crack opening stress 
calculated with equation (32). The predictions, in 
terms of the ratio of predicted Np to test Nt cycles, 
are given in table IV. Three methods were used to 
calculate the stress intensity factor: AIBE, AIBE 
with the hole correction factor, and the solution for 
cracks propagating from a single hole in an infinite 
body. The AIBE method accounts for the boundary 
and crack interaction effects. The hole correction 
factor accounts for the effect of the hole without 
explicitly modeling each hole. The stress intensity 
factor solution for a hole neglects all crack interaction 
effects. 

Table IV. Summary of MSD Fatigue Crack Growth 
Life Predictions 


Test 

Test type 

Analysis method 

Hole in 
infinite 
plate 
solution, 
N p /Nt 

AIBE 

only, 

N p /Nt 

AIBE with 
hole 

correction 

factor, 

Np/Nt 

A6-01 

Nearly equal MSD 

1.66 

1.30 

1.20 

A6-02 

Unequal MSD 

1.46 

0.93 

0.93 

A6-03 

Nearly equal MSD 

1.76 

1.24 

1.14 

A6-04 

Unequal MSD 

2.10 

1.13 

1.13 

A6-05 

Unequal MSD 

1.47 

0.94 

0.94 

Average 

1.69 

1.11 

1.07 


The crack growth behavior for the two MSD tests 
with nearly equal initial crack lengths (no saw cut, 
A6-01 and A6-03) are shown in figures 19 and 20. 
The crack growth behavior of each of the 20 crack tips 
was predicted and the average crack length plotted 
as a function of the applied cycles. The fatigue crack 
growth predictions using the three stress intensity 
factor calculation methods are shown in both figures. 
Neglecting the crack interaction effects resulted in 
a prediction of the crack growth lives that was as 
much as 75 percent greater than observed. The AIBE 
method overpredicted the fatigue crack growth lives 
by as much as 30 percent, primarily because the 
holes will elevate the stress intensity factors when 
the cracks are small. The AIBE method with the hole 
correction factor gave the best predictions, predicting 
lives that were only 20 percent greater than the 
experimental results. 



Figure 19. Predicted and actual fatigue crack growth behavior 
for MSD cracking test A6-01 (nearly equal MSD cracking). 



Figure 20. Predicted and actual fatigue crack growth behavior 
for MSD cracking test A6-03 (nearly equal MSD cracking) . 



Figure 21. Predicted and actual fatigue crack growth behavior 
for MSD cracking test A6-02 (unequal MSD cracking) . 
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AIBE with 



Figure 22. Predicted and actual fatigue crack growth behavior 
for MSB cracking test A6-04 (unequal MSD cracking). 



Figure 23. Predicted and actual fatigue crack growth behavior 
for MSD cracking test A6-05 (unequal MSD cracking). 


The predictions of the crack growth behavior for 
the three MSD tests containing saw cuts to simulate 
crack linkup (unequal initial crack lengths, A6-02, 
A8-04, and A6-05) are shown in figures 21-23. The 
fatigue crack growth predictions made from the three 
methods are shown in comparison with the test re- 
sults. The fatigue crack growth rates for the unequal 
MSD cracking tests were governed by the growth of 
the longest crack. The crack length was such that 
the hole had only a small effect on the stress inten- 


sity factor; thus, the AIBE method and the AIBE 
method with the hole correction factor gave about 
the same life. Neglecting the crack interaction effects 
overpredicted the lives by as much as 100 percent, 
primarily because of the crack interaction and linkup. 
The ratio of predicted to test life ( N p /Nt ) is given in 
table IV for each of the tests. 

Conclusions 

The ability to evaluate the structural integrity of 
aircraft structures with multiple-site damage (MSD) 
is of great concern to the commercial transport air- 
craft industry. An efficient technique was developed 
to calculate stress intensity factors for multiple in- 
teracting cracks. The alternating indirect boundary 
element (AIBE) accurately models multiple interact- 
ing cracks and approximates the local structural de- 
tails with a correction factor. This approach was 
verified through comparison of stress intensity fac- 
tors calculated from accepted solutions, finite ele- 
ment methods, and boundary force analyses. The 
AIBE technique, in conjunction with a correction 
factor for cracks propagating from holes, was used 
to predict the fatigue crack growth life for multiple 
cracks propagating from a line of holes. The predic- 
tions were compared with experimental observations 
and the following conclusions were made: 

• The AIBE analysis is an accurate method of 
calculating the stress intensity factors for multiple 
interacting cracks. 

• The use of the AIBE analysis in conjunction with 
a local correction factor will result in accurate 
stress intensity factors for multiple interacting 
cracks propagating from local structural details 
such as holes. 

• The fatigue crack growth life predicted from the 
AIBE analysis with a hole correction factor was 
within 20 percent of the experimentally observed 
life of tests conducted on flat panels with multiple 
cracks propagating from open holes. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
July 9, 1992 
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